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The factorisation method for Schrodinger operators with magnetic fields on a two- 
, dimensional surface M 2 with non-trivial metric is investigated. This leads to the new 

Q-C integrable examples of such operators and brings a new look at some classical problems 

pH , such as Dirac magnetic monopole and Landau problem. The global geometric aspects 

and related spectral properties of the operators from the factorisation chains are dis- 
cussed in details. We also consider the Laplace transformations on a curved surface 
and extend the class of Schrodinger operators with two integrable levels introduced in 
the flat case by S.P.Novikov and one of the authors. 



Abstract 



X. 

1 Introduction 



In spite of the fact that some important examples of integrable Schrodinger equations with 
magnetic fields were known since 1930-th (Landau problem, Dirac magnetic monopole), 
the general problem of integrability for such equations is still far from being understood. 

In two dimensions probably the first important step in this direction has been done 
by Dubrovin, Krichever and Novikov in 1976, who introduced a very important class of 
Schrodinger operators with magnetic fields integrable ("finite-gap") on one energy level 
[Qj. The coefficients of the corresponding operators are periodic (or quasiperiodic) so that 
the total magnetic flux is zero. 

The case of periodic magnetic fields with non-zero flux has been considered by Dubrovin 
and Novikov || following Aharonov-Casher observation || that the Pauli operators for 
spin 1/2 particles in magnetic field are related to the factorisable Schrodinger operators. 
This allowed to describe explicitly the ground states of the corresponding operators (see 
MM for the details). 
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In H Novikov and one of the authors found a class of operators with magnetic fields 
which are integrable on two different energy levels including the ground state. In section 
5 of this paper we present some generalisations of this result, but our main goal here is to 
investigate what the factorisation method can give for the theory of Schrodinger operators 
L with magnetic fields on a curved two-dimensional surface M 2 . 

Although some ideas have been developed already in the 19-th century by Darboux, 
Moutard et al (see e.g. [jjj), it was probably Schrodinger who first used the factorisation 
method in quantum mechanics Q (see also Q). The idea of this method is very simple: 
if a given operator L can be factorised as 

L = D X D 2 

then the new operator 

L = D 2 D 1 

has the same spectrum as L (provided the operators are good enough) except possibly 
A = and if one knows the eigenfunctions ip of L then the formula 

tp = D 2 ip 

gives the eigenfunctions of L. One can obviously include a shift into the factorisation 
scheme: 

L = D X D 2 + C ^L = D 2 D\ + c 

where c is a constant. For one-dimensional Schrodinger operator it is always possible 
to continue this procedure infinitely and construct a chain of operators related by this 
transformation usually called the factorisation chain (or dressing chain). In the case 
when this chain is periodic the spectrum and eigenfunctions of all these operators can 
be described explicitly (see ||Io|j ), The classical example is harmonic oscillator when the 
corresponding L = L + const. 

In two dimensions the factorisation of the general Schrodinger operator L with mag- 
netic field on a curved surface M 2 is possible only in the special case when the potential is 
equal up to a sign to the magnetic field (see Theorem 1 below). A simple but important 
calculation shows that the transformation 

L = D*D -> L = DD* 

changes the magnetic field B by the Gaussian curvature K of M 2 : 

B = B + K. 

In contrast to the one-dimensional case in two dimensions one can not in general continue 
the factorisation procedure. We give a complete classification of all possible factorisation 
chains on a curved surface M 2 (Theorem 2) which in particular says that infinite factori- 
sation chain exists only on the surfaces with constant Gaussian curvature. The magnetic 
field and the potential of the corresponding operators must also be constant. 

In section 3 we discuss the global geometric aspects of the factorisation chains. The 
Index theorem and the classical Gauss-Bonnet formula play here a crucial role. They 
explain a big difference between positive and negative curvature cases. In the positive 
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constant curvature case when M 2 = S 2 is the standard round sphere we have the Dirac 



magnetic monopole [11], [12], [13]. We show how factorisation method leads to the com- 
plete description of the spectrum of the corresponding Schrodinger operator L in the same 
way as Schrodinger did for the harmonic oscillator in ||. For the flat torus we have the 
standard Landau problem [O]. We also discuss what the factorisation method gives for the 
analogue of the Landau problem on a surface of constant negative curvature with genus 
more than 1. 

In section 4 we present some new examples of integrable Schrodinger operators with 
magnetic fields related to the two-term factorisation chains on the surfaces with non- 
constant curvature. The main observation here is that if the Laplace-Beltrami operator 
on M 2 is integrable then the same is true for the operator with the additional magnetic 
field B = ±K and the potential U = K, K being the Gaussian curvature. 

In the last section we consider the quasi-cyclic chains of the Schrodinger operators with 
magnetic fields on a curved surface related by Laplace transformations generalising the 
constructions from the paper ||. This leads to a class of the operators with two known 
energy levels, one of which is the ground state. 

2 Factorisation method on curved surfaces: local theory 

Consider an oriented analytic surface M 2 with a Riemannian metric ds 2 . It is known that 
M 2 has a complex structure such that the metric is conformal. In any complex chart 
z = x + iy, z = x — iy the metric has the form 

, 9 dzdz dx 2 + dy 2 
ds 



h 2 (z,z) h 2 (z,z) 
The Laplace-Beltrami operator can be defined locally as 

A h = Ah 2 dd = h 2 {d 2 x + d 2 ) 

where d = d z = \(d x —id y ), 8 = ds = \(d x + id y ). To introduce magnetic field one should 
replace the usual derivatives by their covariant counterparts: 

V x = d x — ia, V y = d y — ib, 

V = d-iA, V = d-iA, (1) 

A = ~{a-ib), A = \{a + ib). 
The corresponding Schrodinger operator 

L A = -h 2 \{8 X - ia) 2 + (d y - ib) 2 ] 

can be rewritten as 

L A = -4 h 2 VV + h 2 H = -4 h 2 VV - h 2 H = -2 h 2 (VV + VV) (2) 

where 

H = d x b -d y a = i [V x ., Vy] = 2 [V, V]. 



3 



Geometrically, we have a complex [/(l)-bundle over M 2 with the connection form 

a = i (adx + bdy) (3) 

and the curvature 

f2 = da = i H dx A dy = — — dz A dz (4) 

(see the next section for further discussion of the geometric aspects). Let's define magnetic 
field B by the relation 

n = iB da (5) 
where da = dx Ady is the area element of the surface. By definition we have 

B = h 2 H (6) 

The most general Schrodinger operator on M 2 has the form 

L = -h 2 \{8 X - ia) 2 + (d y - ib) 2 } + U, (7) 

where the potential U(x,y) is a real function on M 2 . 

Consider the following factorisation problem for (|7]): when L can be represented locally 

as 

L = { ai d + a Q ){a\B + a* ) (8) 

(a-factorisation) or 

L = (/3 1 a + /3 )(/3 1 *<9 + /? *) (9) 

(/3-factorisation) for some functions a, a* , (3, (3* ? In the Euclidean case such factorisations 
appeared in the theory of Pauli operators for spin 1/2 particles (see ||, 0]). On a curved 
surface the situation is pretty similar. 

Theorem 1 The a-factorisation for the Schrodinger operator ffi) exists iff U = —B. 
Similarly, the necessary and sufficient condition for the (3 -factorisation is U = B. 

The sufficiency readily follows from formulae (||): 
ifU = -B then L = —Ah 2 VV; 
if U = B then L = —Ah 2 VV. 

The factorisation is not unique: if L = D1D2 is any factorisation, then L = (D\f^ 1 )(f D2) 
is another one for an arbitrary function /. This actually gives all such factorisations. 

Having a factorised operator L = D1D2 one can consider the new operator 

L = D 2 Dx. 

Notice that the change of factorisation D% — > Dif^ 1 , D2 — > fD 2 corresponds to the gauge 
transformation of L: 

L - fD 2 D 1 f- 1 = fLf- 1 . 

The magnetic field and the potential do not depend on the gauge and are defined correctly. 
Let's compute them for a-factorised L. We can assume that D 2 = V, D% = —4 h 2 V so 
that 

L = D 2 D X = -4 V (/i 2 V) = -4 h? VV - 2 h- z h V = 

(10) 

-4 h 2 (V + 2h- z h- x ) V = -4 h 2 VV, 
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where V = V + 2 hgfT 1 , V = V. 

Remark. The nonsymmetry between V and V can be easily corrected by a suitable 
gauge transformation: V — > /iV/i _1 = V — h z h~ l , V — > /iV/i _1 = 'V + hgh . In the future 
we will not worry about such a nonsymmetry provided the corresponding magnetic field 
is real. 

The new magnetic field is 

B = 2 h 2 [V, V] = 2 h 2 [V, V + 2 ^ h' 1 } = 

B + Ah 2 (h- z h~ x ) z = B + 4 ft 2 (In = £ + /i 2 A In h = B + K, 

where K is the Gaussian curvature of the surface. Here we have used the standard formula 
for K in the conformal coordinates 

K = h 2 Alnh, (11) 



(see e.g. [|15|). Thus the new magnetic field is 

B = B + K. (12) 

According to the Theorem 1 the new potential is 

U = B = B + K. (13) 

Similarly, for /3-factorisation the new magnetic field and the new potential are 

B = B — K, U = K - B. 

Let's introduce the notation (B, U) for the gauge class of Schrodinger operators (|7|) with 
given magnetic field B and potential U. Then we have the following two-term factorisation 
chains: 

{B,-B) A (B + K,B + K) (14) 

and 

(B,B) A (B - K,K - B). (15) 
Notice that aof3 = (3oa = Id: 

(B, -B) ^(B + K,B + K)^ (B, -B), 

(B, B) A (B - K,K - B) A (B, B). 

The question is: can we continue the a-chain allowing a shift by a constant? The answer 
is simple: in order to have a-factorisation for an operator from the class (B + K,B + K — 
2c), c = const we should require B + K — 2c = — (B + K) or B + K = c. In this case we 
have the following 3-term chains with an arbitrary constant c: 

(c - K,K - c) A (c, c) = (c, -c) + 2c A 

(16) 

A (c + if, c + K) + 2c = (c + A", 3c + K) 
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and similarly 

(c + K, c + K) ■£ (c, -c) = (c, c)-2c^ 

(17) 



^ (c - K, K - c) - 2c = (c - K, K - 3c). 

Notice that the operators with constant magnetic field and constant potential can be 
considered as the natural analogues of the Landau operators on a curved surface M 2 . 
One more step in the factorisation procedure is possible only if the Gaussian curvature is 
constant: K = Kq. In that case we actually can perform infinitely many steps 



(c-K , K -c)^(c,c)^(c + K , 3c + K )^. 

A ( c + mK , (2m + l)c + m 2 K ) A ... 
for any m G Z + and similarly 

(c + K , c + K )^(c,-c)^(c-K , K -3c)^ 

A (c - mifo, ~(2m + l)c + m 2 if ) .... 
Summarising all this we have 



(18) 



(19) 



Theorem 2 All possible factorisation chains on a surface M 2 with Gaussian curvature 
K are given by 

1) two-term chains (14); (15) with an arbitrary magnetic field B ; 

2) three-term chains (It), (17) of the operators on M 2 with constant magnetic fields. 
If a surface M 2 has a constant Gaussian curvature then we have also 

3) infinite chains tU]l), (Ji^) with constant magnetic fields. 



3 Global geometry and spectral properties of L 

Let's assume now that M 2 is a closed surface of genus g with given Riemannian metric, 
£ is a complex £/(l)-bundle over M 2 , a is a connection on £, Q = iB da is its curvature 
form. Here as above da is the area element of the surface determined by the metric and 
B is a function on M 2 called magnetic field. 

Remark. We should mention that for a given magnetic field B the corresponding 
connection a is defined uniquely modulo natural gauge transformations only when M 2 is a 
topological sphere. For a surface of genus g one can always add to a a closed 1-form without 
changing B. Modulo exact forms corresponding to the gauge transformations these forms 
form the first cohomology group H (M 2 ,R) ~ R 2g . These additional 2g parameters are 
called the Aharonov-Bohm fluxes (see e.g. [jig])- All the spectra we consider in general 
depend not only on the magnetic field B but also on the choice of the connection a and 
therefore on these parameters. When g = 1 (i.e. when M 2 is a torus) this corresponds to 
the choice of 2 Bloch quasi-momenta. 

Notice that the total magnetic flux is an integer multiple of 2ir: 

— ( B da = b, b G Z. (20) 
2tt Jm 2 
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The integer b is actually the first Chern class of £ : 



b = ci(£) G H 2 (M 2 ,Z) ~ Z (21) 

(see e.g. |15[| ). Having all this plus a potential [/ which is a function on M 2 one can 
define the Schrodinger operator L ([?]) as explained in the previous section. This operator 
is acting on the sections S(£) of the bundle £ and is self-adjoint with respect to the natural 
Hermitian structure with the norm 

IMI 2 =/ \i>\ 2 do. 

Jm 2 

Let's analyse now what happens when we apply the factorisation procedure. Consider 



first the two-term chain (14): 

(B,-B) A (B + K,B + K) 

or in the local coordinates 

L = D*D -> DD* = L 

where D = 2V, D* = -2/i 2 V. Here we have used the complex structure on M 2 uniquely 
determined by the Riemannian metric, which exists according to the classical results (see 
e.g.®)- 

Notice first of all that the operator D maps the sections S(£) of £ into the sections S(rj) 
of the bundle r] = £ (g) T 0,1 (M 2 ) which are the antiholomorphic 1-forms on M 2 with the 
values in £ of the form ipdz, tp G S(£). This space also has the natural Hermitian structure 
induced by the Hermitian structure on £ and the metric on M 2 . It is easy to check that 
D* = —2h 2 V is indeed operator adjoint to D = 2V with respect to these structures, so 
the operator L = D*D is non-negative and KerL = KerD. 

The new operator L = DD* is acting on the sections S(rj) of the bundle r\ and in the 
local coordinates has the form ([l0[): 

L = -4 V (/i 2 V) = -4 h 2 (V + 2h z h~ 1 ) V. 

The new covariant derivatives V = V + 2 h^h^ 1 , V = V correspond to the natural connec- 
tion on 7] induced by the connection a on £ and the natural Hermitian connection on the 
antiholomorphic cotangent bundle T 0,1 (M 2 ) (see [17|). The curvature of this connection 
is 

Cl = i(B + K)da, 

where K is the Gaussian curvature of the metric. Notice that the total flux of the new 
magnetic field B = B + K (or, equivalently, the first Chern class of the bundle rf) is 

ci(rj)=b=7rf B da = b + X = b+(2-2g) (22) 

2-7T J M 2 

where \ = 2 — 2g is the Euler characteristics of M 2 because of the Gauss-Bonnet formula 

K da = 2tt X . (23) 

M 2 
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The operator D is an elliptic operator from S(£) into S(rj). Its index can be determined 
by the Index theorem and is given by the Riemann-Roch formula 



ind D = dim Ker D - dim Ker D* = b - g + 1. (24) 
When magnetic field is large enough, more precisely when b > 2g — 2, we have 

dim Ker D* = 0, dim Ker D = b - g + 1. (25) 
In that case the ground state of the operator L is degenerate: 

dim Ker L = dim Ker D = b — g + 1 
while L is positive operator: 

dim Ker L = dim Ker D* = 0. 

The rest of the spectrum (which is discrete according to general theory) is the same 
for L and L: the intertwining operators D and D* establish the isomorphism of the 
corresponding eigenspaces. Notice that if b = g — 1 then, according to the Riemann-Roch 
formula (p3|), dim Ker L = dim Ker L so that L and L are isospectral everywhere ( not 
just for A > 0). Now for the general factorisation chain one should only take into account 
the additional shift of the spectrum. 



For the /3-factorisation chains (15) the analysis is similar, one should simply replace 
the holomorphic structure by the antiholomorphic one. 
Let's consider now some examples. 

Example 1. Dirac magnetic monopole on a sphere. 

Let M 2 be a sphere S 2 C R 3 with the standard metric of constant Gaussian curvature 
K. The Hamiltonian of the Dirac monopole with a charge q £ Z in our notations is a 
Schrodinger operator H q from the gauge class (B,0) where B is a constant satisfying the 
quantisation relation 

_L f Bda=^- [ da = 2B = q, 

27T J52 2-7T Js2 

i.e. B = q/2 must be integer or half-integer. It is acting on the sections of the U(l)-bundle 
with the first Chern class q. As we have seen above such operator can be included in the 
infinite factorisation chain 

(B, 0) = {B, -B) + B ^(B + K, B + K) + B = (B + K, 2B + K) = 



{B + K, —{B + K)) + 3B + 2K -» (B + 2K, B + 2K) +3B + 2K = (26) 

(B + 2K, AB + AK)^ ... ->(B + mK, 2mB + m 2 K) -> 

m € Z+. If B = q/2 > then the operator (B, —B) is positive and according to index 
theorem ( |25| ) 

dim Ker (B, -B) =q-g + l=q+l 

since g = and q > 2^ — 1 = —1. Therefore the ground state of H q has the energy 
A = B = q/2 and the corresponding eigenspace has dimension q + 1 = 2B + 1. The second 
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operator (B + K, 2B + K) in the chain has the same spectrum as (B, 0) except the ground 
state. By the same reasons its ground state has energy A = 3B + 2K which is degenerate: 



dim Ker [(B + K, 2B + K) - (3B + 2K)\ = dim Ker [B + K, -(B + K)] = 

2B + 2 + l = 2B + 3. 

Thus the second eigenvalue of H q is A2 = 3B + 2K with degeneracy 2B + 3. On the m-th 
step we'll have the operator (B + mK, 2mB + m 2 K) = (B + mK, —B — mK) + (2m + 
i)B + (m 2 + m)K which leads to the eigenvalue 

A m = (2m + 1)B + m(m + l)K 

with degeneracy 2B+2m+l. Thus we arrived at the well-known result about the spectrum 
of the Dirac monopole (see e.g. [fiT 



Spec H q = {A n = (2m + 1)- + m(m + 1)K with degeneracy g + 2m + 1}, 
m = 0,1,2,... 

We can make all this explicit (including the calculation of the corresponding eigen- 
functions) using the stereographic coordinate z = x + iy. Assume for simplicity that the 
radius of the sphere R = 1, so that the Gaussian curvature K = 1/R 2 = 1. The metric 
has the form ^ 

ds 2 = —77 dzdz. 

(l + zzy 

The corresponding chain of operators in a suitable gauge has the form 

L N = -(1 + zz) 2 dd - Nz(l + zz) d + Nz(l + zz) 8 + N 2 (l + zz) = 

D* N D N + N(N + 1) = D N ^D* N ^ + N(N - 1), (27) 

D N = (1 + zz) d + Nz; D* N = — (1 + zz) d + (N + l)z. 

Here = H2N + N 2 , where H2N is the Dirac monopole operator with the charge q = 2N, 
B = N. We obviously have the intertwining relations 

Ln+\Dn = DnLn, LnD* n = D* n Ln + i. (28) 

To find the ground state of : Lj^i/j = N(N + l)ip one should solve the equation 
D N i) = 0: 

- - Nz 

(1 + zz) dip + Nz ib = 0, or d hi ib = 

1 + zz 

The solutions are easy to find: 



ib 



(1 + zz) N 

where / is any holomorphic function of z. Because of the condition 



1 2 



, , \ib\ dzdz 
da = I I —— — — ^ < ex) 

R 2 (1 + zzy 
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the function f(z) must be a polynomial of order < 2N. This gives us the space of dimension 
2N + 1. Applying to this space the "lowering" operators Z)^ r _ 1 , -D^_ 2 , ... we will construct 
all eigenfunctions of the operators -Ljv-i, £jv-2> ••• with the eigenvalue A = N(N + 1). This 
gives the following description of the eigenfunctions of the Dirac monopole operator H q . 

Theorem 3 The eigenfunctions of the Dirac magnetic monopole operator with the charge 
q on a unit sphere corresponding to the eigenvalue 

A = (2m + 1)~ + m(m + 1), m G Z + (29) 
form the space of dimension q + 2m + 1 which can be described as 

t/j = D* N _ m ... Dj f _ 2 Dj f _ 1 ^ zz) N 
where N = m + % and f(z) is an arbitrary polynomial of degree deg/ < 2N = 2m + q. 

One can check that the formulas for tfi we have given actually determine the smooth 
sections of the corresponding line bundles over the sphere S 2 . 

The eigenfunctions of the Dirac magnetic monopole are known as monopole harmonics 
and have been investigated by Wu and Yang [^], who were probably the first to identify 
them explicitly as the sections. Our derivation is different and more close to the one of the 
paper |l9[| by D'Hoker and Vinet who discovered the supersymmetry of the corresponding 
Pauli equation in the presence of a magnetic monopole (see also |2(| ) . 

Remark. If q = 2N is even integer then H^n = Ln — N 2 can be intertwined with the 
shifted standard Laplace-Beltrami operator Lq = — A g on the sphere S" 2 : 

H 2N D = D(L - N 2 ), D = D N _! ... D X D Q , 

so that one can use the well-known eigenfunctions of Lq (spherical harmonics) to construct 
the eigenfunctions of H q . Notice that we have given an alternative description of spherical 
harmonics using the factorisation chain. The explicit form of the intertwining operator D 
is 

D = t, 9(1 + zz) 2 B(l + zz) 2 5 ... 8(1 + zz) 2 8 

(1 + zz)" z 

which is of the order TV in 8. 

It is interesting to compare the Dirac monopole problem on a sphere with the Landau 
problem on the surfaces of genus g > 2 with the constant negative Gaussian curvature K < 
which we analyse below. The corresponding classical problems behave very differently 
so one should expect the same for the quantum problems as well. 

For the investigation of the intermediate case of flat torus which is the classical Landau 
problem we refer to [14] and ||. In that case K = and our formula (p6|) leads to the 



usual Landau spectrum 

A m = {(2m + 1)B, m = 0,l,...}, 

where B is assumed positive and quantised. The corresponding eigenfunctions can be 
expressed in terms of the classical elliptic <r-functions (see ]|], ||] for the details). 
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Example 2. Landau problem on a surface of genus g > 2. 

Let M 2 be any analytic surface with the metric ds 2 of constant curvature K = — 1, 
Such surface M 2 can be considered as a quotient of the Lobachevsky plane C by an infinite 
discrete group G: 

M 2 = C/G. 

C can be realised as an open disc \z\ < 1 with the metric 

ds 2 = t — tt dzdz. 

(1 — zz) A 

The Schrodinger operator on C with constant magnetic field B (Landau operator) can be 
written as = -(1 - zz) 2 (VV + VV), where 

V = 9- J B-^— , V = 8 + B-^—, 

1 — zz 1 — zz 

or, explicitly, 

L% = -(1 - zz) 2 dd + Bz(l - zz) 8- Bz(l - zz) d + B 2 zz. (31) 

In order to define the corresponding operator Lb acting on the sections of some U(l)- 
bundle over M 2 one needs the quantisation condition 

— / B da = —B [ da = (2o — 2)B £ Z (32) 
2vr i M 2 2vr A/ a 

to be satisfied. Here we have again used the Gauss-Bonnet formula 

— / K da = — - / da = x = 2 - 2g. 

27T Jm 2 2-7T Jm 2 

Notice that to define Lb one also needs to choose a connection a which depends on 2g 
Aharonov-Bohm fluxes (see the remark at the beginning of this section). 

Let's see what the factorisation chain ([0]) gives us for the calculation of the spectrum 
oi Lb- We have the chain (^) with K = — 1: 

(B,0) -> (B,-B) + B -> (S - 1,5 - 1) + B = (B - l,2B - 1) -> ... 
— > (B — m, 2mi? — m 2 ) = (B — m,m — B) + (2m + l)i? — (m 2 + m). 



Let's assume that B > 0. By Riemann-Roch formula (|24j) the index of the operator 
(B — m : m — B) is (2g — 2)(B — m) — g + 1 = (2p — 2) (S — m — ^), so it is positive if 

m < B . 

2 

Using the same arguments as in the previous example we can claim that the first [B— h] + 1 
eigenvalues of the operator Lb have the form 

X m = {2m + l)B -m(m + l), m = 0, 1, 2, [B - -]. (33) 
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Let's remind that B has the form k/(2g — 2) for some positive integer k. Moreover, we 
can say that if (2g — 2)(B — m) > 2g — 2, i.e. if 

m < B — 1 

then the corresponding eigenspace has the dimension 

dim Ker (L B - X m ) = (2g - 2){B - m - -). (34) 

Notice that the spectrum (|33| ) depends only on the magnetic field B but not on the 
Aharonov-Bohm fluxes. It is related to the discrete part of the spectrum of the Landau 
operator on the whole Lobachevsky plane (see e.g.pl|1). It would be interesting to find 
an effective representation for the corresponding eigenfunctions on a surface M 2 given 
explicitly as an algebraic curve in C 2 . 

From these considerations we have nothing to say about the rest of the spectrum 
oi Lb (which in fact depends both on Aharonov-Bohm fluxes and on (3g — 3) complex 
parameters (moduli) determining the conformal structure on M 2 ). In particular, we can 
not say anything about the spectrum of the pure Laplace-Beltrami operator Lq on M 2 . 
This is a reflection of the fact that the corresponding classical geodesies problem on M 2 
is non-integrable. 



4 New examples of the integrable quantum problems with 
magnetic field 

Let's consider any surface M 2 with integrable quantum geodesic problem 

Lip = Xip, 

where L = —A^ and A/j is the Laplace-Beltrami operator on M 2 . By definition this means 
that there exists another differential operator F which commutes with L: [F, L] = and 
has an independent highest symbol. We have shown (see ( |l6|) above) that L can be 
factorised L = D*D and the new operatorL = DD* has the magnetic field B = ±K (K 
is the Gaussian curvature) and the potential U = K. We claim that this new operator is 
also integrable. Indeed, consider the differential operator 

F = DFD* , (35) 

then 

LF = DD*DFD* = DLFD* = DFLD* = DFD*DD* = FL, 

i.e. F commutes with L. It is easy to check that if the highest symbol of F is independent 
of the highest symbol of L then the same is true for F and L. The intertwining relations 

LD = DL and D*L = LD* 

establish isomorphism between the two spectral problems 

Lip = Xip and Lip = Xip 
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for any A / : ip = Dip, ip = D*ip. For A = the situation is described by the 
Riemann-Roch formula: 

dim Ker L — dim Ker L = 1 — g (36) 
In particular, for the sphere S 2 we have g = 0, 

dim Ker L = 1, dim Ker L = 0. 

This means that the ground state of L has the energy which is equal to the minimal 
positive eigenvalue of L: if Spec L = {Ao = 0, Ai, A2, •••}, then Spec L = {Ai, A2, •••}• 
For the torus T 2 we have g = 1 and 

dim Ker L = dim Ker L = 1 



so that 

Spec L = Spec L. 

This is true for all the values of quasi-momenta i.e. for the whole Bloch spectrum. Notice 
for the torus the total magnetic flux of L is zero. 

We do not know any integrable Laplace-Beltrami operators on a surface of genus g > 2. 
For the corresponding classical problem about geodesies on M 2 there exists a rigorous proof 
that there are no such metrics (see |16||), 

Theorem 4 Let M 2 be any surface such that the corresponding quantum geodesic problem 
—Afrip = Xip is integrable. Then the Schrddinger operator L on M 2 with magnetic field 
B = ±K and potential U = K, K is the Gaussian curvature of M 2 , is integrable too and 
has the same spectrum as L = — with the only possible exception at A = 0. 

Example. Consider an ellipsoid M 2 

x 2 y 2 z 2 
V - — I = 1 

a 2 b 2 c 2 

with the metric induced from R 3 . The geodesic problem on M 2 has been solved by Jacobi 
who showed that it can be integrated by separation of variables. The same is true for the 
corresponding quantum problem. Gaussian curvature of M 2 has the form 

so we can claim that the Schrddinger operator L on M 2 with the magnetic field 
and the potential 
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is integrable and Spec L = Spec L U {0}. Notice that the order of the additional quantum 
integral F = DFD* is 4, since the order of F is known to be 2. 

Remark. We do not claim that that the minimal order of the additional integral is 
4. On the contrary, in this case one can show that there exists an additional integral 
of order 2. It would be interesting to investigate the corresponding classical mechanical 
problem of motion on the ellipsoid in this special magnetic field. We conjecture that it is 
nonintegrable. 

This example can be generalised in the following way. Consider any surface M 2 with 
the Liouville metric 



ds 2 = gndu 2 + g 22 dv 2 , gn = — — , g 22 

J g 

where f(u) and g{v) are arbitrary functions of the specified arguments. Its Gaussian 
curvature is expressed by the formula 

K f-g f' + g' 



u — V v — u 



2(u-v) 3 A(u-v) 2 ' 
The corresponding Laplace-Beltrami operator 



Q n I Q 22 



commutes with the second-order operator 



g 11 Q , .. rrr^ Q g 22 



F = v^g^g 22 d u -j== d u + u^Jg^g 22 d v -j== d v , 

so that variables can be separated both in the classical and quantum cases. So we can 
claim that the Schrodinger operator with magnetic field B = ±K and the potential U = K 
is integrable on any Liouville surface. 

As a degenerate case of this construction one can get the surfaces of revolution in R? 
with the metric 

ds 2 = p{z) 2 d^ 2 + (1 + p'{z) 2 ) dz 2 . 
In that case the Gaussian curvature is given by the formula 

K 



P(l + W) 2 ' 

We would like to mention that there exist the intertwining relations between two 
Schrodinger operators which are not related to any factorisation of the operators (see |22j 
for some examples). We will present here another example of this type modifying the 
previous analysis of the Dirac monopole. Recall that the operator of Dirac monopole 
on the unit sphere given by (p7|) 

L N = -(1 + zzf dd - Nz(l + zz) d + Nz(l + zz) 8 + iV 2 (l + zz) 

satisfies the interwining relations 

Ln+\F>n = DnLn 
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where 

D N = {1 + zz) d + Nz. 
It can be readily verified that the modified operators Ln given by 

In = L N - i P 2 - (N + 1) P, L N+1 = L N+l -^P 2 -NP 

satisfy the same intertwining relations 

Ln+\Dn = DnLn 
where the modified intertwining operator is of the form 

D N = D N + r 



1 + zz 

Here ip is a quadratic polynomial in z: 

ip = q+pz-qz 2 , p G R, q,q~eC, 
and P is given by the formula 

, 2zv? 1 - ^ + qz 

P = tp z = p z — 2 — . 

1 + zz 1 + zz 1 + zz 

Geometrically, P represents the restriction to the sphere S 2 of an arbitrary linear function 
from the ambient space i? 3 . Without any loss of generality we may assume q = q = (by 
appropriately choosing the axis of stereographic projection), so P = p j^ff- 

In particular for iV = lwe have the intertwining relation between the Dirac monopole 
operator with charge 2 in the potential U = — \ P 2 and the usual Laplace-Beltrami 
operator with the additional potential V = — \ P 2 — P. 

One can show that this construction actually gives all the potential deformations of 
the intertwining relations (|2"8|). 

5 Laplace transformations on a curved surface and quasi- 
cyclic chains 

Let L be any Schrodinger operator with magnetic field B and potential U: 

L = -4 h 2 VV + (U + B) = -4 h 2 VV + (U - B). 

In general none of these forms is pure factorisable but on the level A = : Lip = we 
still can do the transformation 

ip = VV' or i\) = Vifj. 

In particular, if 

Lip = [-4 h 2 VV + (U — B)] ip = 

and if) = \7ip we have 

U-B , 
V^ = ^^ 
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and therefore 

U-B ~ 
VW = ^^ 

where V satisfies the relation 

~U -B U -B „ 

V —TK = V 

Ah 2 Ah 2 

implying 

V = V - ln(C7 - B) z + (In h 2 ) z . (37) 
Thus t/j = satisfies the new Schrodinger equation Lip = 0, where 

L = -A h 2 VV + (U - B) = -A h 2 VV + (U + B). (38) 

The new magnetic field is 

B = 2 h 2 [V, V] = B + 2 h 2 Qn(U-B)) zS -4: h 2 Qnh) zg = B + - A h (ln(U - B)) - K (39) 

where Ah = Ah 2 dd is the Laplace-Beltrami operator on M 2 , K is the Gaussian curvature. 
The new potential is 

JJ = U-B-B = U-[2B + ^ A h (ln(U - B)) - K\. (40) 

The formulae ([??], 38.39,^) define correctly the Laplace transformation for the Schrodinger 
operators on a curved surface. The only difference with the standard formulae in the flat 
case (see e.g. ||) is the additional Gaussian curvature term. 
Following H let's call the chain of Laplace transformations 

B k+ i = B k + \ A h (\n(U k - B k )) - K, 

(41) 

£4+1 = U k - B k - B k+1 = U k - [2B k + i A h (\n(U k - B k )) - K], 

k = 0, N, quasi-cyclic if both the initial and final Schrodinger operators are factorisable, 
possibly at the different energy levels: 

U + B = 0, U N + B N = U N _! - B N ^ = -c (42) 

where c > is a constant. 

As well as in the flat case [|| the last operator Ln has two "integrable" levels: L^ipQ = 
and L]yip c = —c tp c . Indeed, i/jq can be found as the result 

ipQ = V7V_iVat_ 2 ...VoV ; 

of the Laplace transformations applied to the solutions of the initial equation L^ip = 0, 
which is equivalent to 

V ^ = 0, (43) 

while ip c are the ground states of the operator Ljy = —4 h 2 Vat Vat — c satisfying the 
equation 

Vtv^ = 0, . (44) 
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Obviously we should assume that solutions of both (|43|) and ( f44| ) do exist, which imposes 
some global assumptions on the magnetic field. Let's discuss these assumptions. 

First of all let us notice that the magnetic charge b = j- J M 2 B da (which should 
be integer because of the quantisation condition) changes under Laplace transformations 
according to the formula 

b = b + 2g -2 



as it follows from (39) and the Gauss-Bonnet formula. After N steps we have 6/v 



b + 2N(g - 1). Also from @ we have 

/ (U k + B k ) da= ( (C/ fc _i - B fc _0 da = [ (U k ^ + B fc _ a ) da - ^b k _ 
Jm 2 Jm 2 Jm 2 

Since Uo + Bq = it follows that 

l - [ (U N + B N ) da = ~ [ cda = -2Nb-N(N-l)(2g-2), 
!7T J m 2 2tt Jm 2 



1- 



2tt 

so we have the relation 

A(M 2 



2vr 



2Nb + N(N-l)(2g-2), (45) 



where ^4(M 2 ) = J M2 da is the area of M 2 . For the sphere we have c = A r 6o — N(N — 1), 
for the torus cA(M 2 ) = AirNbo, for a surface of genus g > 1 the relation c(2g — 2) = 
2A% + A r (A r -l)(23-2). 

This determines the constant c in the quasi-cyclic chain if we know the magnetic charge 
bo of the first operator and therefore imposes a quantisation condition on c since bo is an 



integer. In terms of this integer &o the sufficient conditions for the equations (43) and ( [44]) 
to have a solution have the form (see the section 3) 

b > g-l,b + 2N(g-l)> g-1. 

For the topological sphere S 2 this is equivalent to the inequality 6o > 2N — 1, for the torus 
this simply means that &o is positive, for a surface of genus g > 1 we have just the first 
inequality bo > g — 1. 



Theorem 5 Let U k ,B k , k = 0, 1, N satisfy the quasi-cyclic chain (41 ) on a curved su- 
face M 2 and let the magnetic charge of the first operator bo satisfy the conditions described 
above. Then the last operator of the chain Ln has two known energy levels independent 
of the Aharonov-Bohm fluxes: the ground state A = — c and A = 0. The corresponding 
eigenf unctions can be found from the solutions of the first order equations and $4% ) 
and for a large bo form the spaces of the dimensions bo + (2N — l)(g — 1) for A = — c and 
^0 — 9 + 1 for A = respectively. 

Example 1. Quasi-cyclic chains of length N = 1: 

Uo + £ = 

U 1 + B 1 = Uo-B Q = -c 
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implies that Uq = — S = — c/2 are constants, so that L is the Landau operator on a 
curved surface M 2 . Notice that 

Si = B + 1 A /l (ln(-c)) -K = B -K, 

Ui = U - 2S + K = K - 3S , 

so that in this case Laplace transformation coincides with one step of the factorisation 
procedure (|l7|): (So, — So) — > (Bq — K, K — 3So), So = §. Thus in this case we claim that 
the Schrodinger operator with the magnetic field S = | — K and the potential U = K — ^ 
has two lowest energy levels known: A = — c and A = provided the quantisation and 
positivity conditions for c are satisfied. 

Example 2. Quasi-cyclic chains with N = 2 : 

U Q + S = & U = -S ; 

Sx = S + iA^lnSo) -K, U l = U Q -B Q -B 1 = -2S - B x ; 

U 2 + B 2 = U 1 -B 1 = -2S - 2Si = -c. 
Thus we have the following relation for the magnetic field So: 

2So + ^A h (lnS )-^=| 

or 

A^lnSo) =c + 2if-4S 
which after introducing (p = In So takes the form 

A h <p = C + 2K-A e^. (46) 

When K = this reduces to the equation from ||. 

As well as in the flat case any solution of ( ]46| ) determines a Schrodinger operator in 
magnetic field with two integrable levels. 

Remark. When c = equation ( 46 ) reduces to 

A h ip = 2K - 4 e v 

which is a natural analogue of the well-known Liouville equation for a curved surface. 
It can be transformed into the standard Liouville equation (p Z z = —e^ by a substitution 
ip = ip + 21n/i. When c / the equation ( f46[) is probably non- integrable already in the 
flat case. 
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